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^ '■ ABSTRACT 

We summarize a general formulation of particle propagation in fluctuat- 
ing media, as applied to the description of neutrino propagation through the 
sun. It contains the familiar MSW effective hamiltonian, plus corrections 
which describe neutrino interactions with fluctuations in the medium. An 
estimate of the size of these corrections for a simple model of solar fluctua- 
tions is made, with the conclusion that they can produce surprisingly large 
effects since they grow with the correlation length of the fluctuation. For 
■ MSW oscillations, the leading effect for small fluctuations is to diminish the 

quality of the resonance, making the suppression of the 7 Be neutrinos an 
experimental probe of fluctuations deep within the sun. Fluctuations can 
also provide an energy- independent suppression factor of ^, away from the 
resonant region, even for small vacuum mixing angles. To be even potentially 
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\ detectable, density fluctuations must be correlated on scales of hundreds of 

kilometres, and have amplitudes of a few percent or larger. 

The next generation of solar-neutrino detectors, including SNO 
and Super Kamiokande,a bring us into a new era, in which neutrino 
solutions to the solar-neutrino problem will be put under detailed ex- 
perimental scrutiny. It therefore behooves theorists to understand in 
detail the corrections to the predictions of standard solar models, with 
and without neutrino oscillations. Since standard treatmentsB of MSW 
oscillations replace the solar medium with an average electron den- 
sity, one such correction consists of how fluctuations in solar properties 
about this average can affect neutrino propagation. We very briefly 
review here a methoda for treating such fluctuations systematically. 
Limitations of space preclude our discussing related approaches^ to 
these issues. 



*Talk presented to Neutrino 96, Helsinki, Finland, June 1996. 
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A great deal is known about how particles like photons and neu- 
trons propagate through materials, and general techniques have been 
developed to describe this propagation. Our approach is based on an 
adaptation of these ideas to neutrino physics. In order to keep approx- 
imations explicit, we start with the general case, and then specialize to 
neutrinos within the sun. In its most general context we consider two 
sectors, A and B, where A describes the degrees of freedom we wish to 
follow - such as neutrino flavour - and B describes the other degrees 
of freedom, which we do not intend to measure - like those of the solar 
medium. (There is a generalization of the approach for which a limited 
measurement to B is also performed.@) We imagine B to be essentially 
unchanged by its interaction with sector A. All predictions for mea- 
surements purely within sector A are described by the reduced density 
matrix, defined by tracing the full density matrix over the degrees of 
freedom in B: p A (t) = Tr B p(t) . 

The time evolution of p A is computed using the hamiltonian H = 
Hq + V, where H = H A + H B describes the separate evolution of sectors 
A and B, and V describes their mutual interaction. There are two 
important timescales which govern this time evolution, and considerable 
simplification arises when they are very different. The first scale, r p , 
is defined as the time beyond which perturbation theory in V fails. 
The second scale, r c , is the correlation time, defined as the time above 
which (5V(t)5V(t')) becomes negligible (if such a time exists). Here the 
average, (■ ■ ■) = Tr s (/>••-), is over sector B, and SV = V — V is the 
deviation of the interaction hamiltonian from the effective interaction, 
V, which describes the average evolution in sector A. In perturbation 
theory V(t) = (V(t)) - i ft dr (5V(t) 5V(t)) + O (V 3 ). 

If t c <C r p , then p A satisfies a master equation which can be de- 
rived in perturbation theory, and then integrated to obtain its large- 
time behaviour, even for t > r p . It is given, to second order in V, by: 

^ = -i[V(t)p A (t)- PA (t)V*(t)] (1) 

+ y* dr Tr B (6V(t) p(t) SV(t) + SV(r) p(t) SV(t)) + O (V 3 ) . 

These ideas may be applied to obtain an equation describing the 
evolution of neutrino flavour within the solar mediumB The result is a 
two-by-two matrix equation: 

?PL = -i[ Vo + V 1 ,p f \ -2G 2 F A[(g e ) 2 Pf + P f (g e ) 2 -2g e p f g e ], (2) 

where pf is the 2x2 neutrino-sector density matrix in flavour space, 
Vq ~ k + + • • •, where k is the neutrino momentum and m is 
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its mass matrix. Vi(t) = y/2 G F (g e n e (t) + g n n n {t)) is the first-order 
effective interaction, with g e = diag(l,0), and g n = diag(— |, — 5) 
being the 2x2 matrices which describe the charged-current (neutral- 
current) neutrino coupling to the electron (neutron) density, n e (n n ). 
Because g n is proportional to the unit matrix, the neutron density drops 
out of eq. (H) . The first term of eq. (@) gives the usual MSW evolution. 
All second-order terms are proportional to the coefficient: 

A{t) = f dr (5n e {t) 5n e {r)), (3) 

which controls the size of fluctuation effects. 

It is instructive to compare the order of magnitude of the various 
terms in eq. (|2]). The flavour-changing part of Vq is m^m/k; while the 
term involving V\ is ~ G F g e {n e ). These two terms are comparable over 
part of the solar interior, and their interplay gives rise to the MSW 
resonant oscillation. The contribution of the remaining term is of order 
G F (5n e 5n e )r c , which is competitive with the others if G F (5n e b~n e )r c ~ 
(n e ). Writing (5n e 5n e ) = e 2 (n e ) 2 , this condition becomes e 2 G F (n e )r c ~ 
1. Since, deep within the solar interior, G F (n e ) < (100 km) -1 , signifi- 
cant fluctuation effects require e 2 r c > 100 km. 

These estimates are borne out by more detailed calculations.Q 
When applied to thermal fluctuations, they reproduce the standard re- 
sult for neutrino scattering from thermally-distributed electrons, giving 
a negligible scattering rate within the sun. We have identified only one 
potential source of fluctuations within the sun whose effects may be 
nonnegligible. These fluctuations arise due to the variation, over space 
and time, of the mean electron density within the sun. The propaga- 
tion of any particular neutrino through such variations is perfectly well 
described by the MSW treatment. In general, however, successive neu- 
trinos do not encounter the same electron density profile. This could 
be because different neutrinos are produced at different places within 
the sun, and so pass through different regions while en route to the 
earth. Alternatively, successive neutrinos produced at the same point 
within the sun could pass through different electron densities because 
the density profile itself generally varies in time - such as when he- 
lioseismic waves are present. As a result, the neutrino flux to which a 
detector is exposed can be thought to have experienced an ensemble 
of density profiles, which must be averaged to obtain the integrated 
neutrino signal as seen by a detector on earth. 

We have computed^ how the MSW resonance is altered due to two 
kinds of ensembles of this type which are of particular interest. These 
ensembles differ in the nature of the correlations which they assume. 
One choice assumes fluctuations which are uncorrelated in position 
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space, i.e. 8n e (x) is uncorrelated with 5n e (y) for sufficiently large |x— y|. 
This kind of local fluctuation is typical of many condensed-matter ap- 
plications, and could arise in the sun from a number of sources. The 
second choice consists of fluctuations which are uncorrelated in momen- 
tum space, with Sn e (k) uncorrelated with 8n e (k'), where k ^ k' label 
an appropriate set of density-oscillation modes. This second kind of 
ensemble is intended as a simple model of how neutrinos interact with 
helioseismic waves, since these have been observed to exist in the sun. 

It is straightforward to generalize the MS W analysis to arrive at 
a variation of Parke's formulaQ: 

P^) = \ + (I ~ Pj) A cos2fl m (t') cos29 m (t), (4) 
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in which Pj = exp [- f (^ff^) (^)] is the usual 'jump' probabil- 
ity as one passes through the resonance point, with h is the scale height 
for the electron density and 9 V , 9 m and 5m 2 = va\ — rn? t are the vacuum 
and matter mixing angles, and the difference between the squares of 
the masses, for the two neutrino states. Fluctuation effects enter eq. @) 
through the damping factor A = exp —2 ) sin 29 m (x) . 

Notice that the fluctuations dominantly act to damp the neutrino 
oscillations. This damping is a reflection of the conversion of the in- 
coming neutrino from a pure to a mixed state as it interacts with the 
fluctuations. The damping ruins the M SW resonance condition, thereby 
reducing the effectiveness of the flavour change as the neutrino passes 
through the resonant point. Our numerical results (presented in the 
figures) show that small fluctuations first become noticeable where the 
resonant flavour change is the strongest. Since 7 Be neutrinos are sup- 
pressed by MSW oscillations in this way, a deviation from the MSW 
prediction for the strength of the 7 Be neutrino flux could be evidence 
for electron density fluctuations deep within the solar interior. 

For sufficiently large times, eq. (H) has the universal prediction: 
P e — > |. This limit is also seen in our numerical integrations, which 
agree well with eq. (^) throughout parameter space. This suggests a 
new solution for the solar neutrino problem: an approximately energy- 
independent suppression of the solar neutrino flux by a factor of 2 due 
to decoherence by solar fluctuations. Such a suppression would look 
much like large-angle vacuum oscillations, even though it could arise 
through fluctuations with small vacuum mixing angles. 

The figures present a summary of the results of a numerical inte- 
gration of eq. (H), using eq. (|j) for density fluctuations which are uncor- 
related in either position and momentum space. As is seen from these 
figures, deviations from MSW predictions require fluctuations whose 
amplitude is, at the very least, a few percent. 
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Figure 2 
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The v e survival probability as a function of the neutrino energy, in 
units of Am 2 . Figure 1 (Figure 2) uses fluctuations which are uncor- 
rected in position (momentum) space. The vacuum mixing angle is 
taken to be 0.1. epsilon denotes the fractional amplitude of the fluc- 
tuation while in Figure 1 £ represents the correlation length. 
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